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Abstract. Simple analytical models, such as the Hernquist model, are very useful tools to investigate the dynamical 
structure of galaxies. Unfortunately, most of the analytical distribution functions are either isotropic or of the 
Osipkov-Merritt type, and hence basically one-dimensional. We present three different families of anisotropic 
distribution functions that self-consistently generate the Hernquist potential-density pair. These families have 
constant, increasing and decreasing anisotropy profiles respectively, and can hence represent a wide variety of 
orbital structures. For all of the models presented, the distribution function and the velocity dispersions can 
be written in terms of elementary functions. These models are ideal tools for a wide range of applications, in 
particular to generate the initial conditions for A^-body or Monte Carlo simulations. 
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1. Introduction 

From a stellar dynamical point of view, the most complete 
description of a stellar system is the distribution func- 
tion F{r,v), which gives the probability density for the 
stars in phase space. In this paper, we will concentrate on 
the problem of constructing anisotropic equilibrium dis- 
tribution functions that self-consistently generate a given 
spherical mass density profile p{r). In the assumption of 
spherical symmetry, the mass density of stellar system can 
easily be derived from the observed surface brightness pro- 
file, at least if we assume that the mass-to-light ratio is 
constant and that dust attenuation is negligible. And as 
the surface brightness of a galaxy (or bulge or cluster) 
is fairly cheap and straightforward to observe, compared 
to other dynamical observables which require expensive 
spectroscopy, the problem we will deal with is relevant 
and important. 

The first step in the construction of self-consistent 
models is the calculation of the gravitational poten- 
tial ^{r), which can immediately be determined through 
Poisson's equation. The second step, the actual construc- 
tion of the distribution function, is less straightforward. 
Basic stellar dynamics theory (see e.g. Binney & Tremaine 
1987) learns that steady-state distribution functions for 
spherical systems can generally be written as a function 
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of binding energy and angular momentum. We hence have 
to determine a distribution function F(£, L), such that the 
zeroth order moment of this distribution function equals 
the density, i.e. we have to solve the integral equation 



p{r) = 



F{£,L) dv 



(1) 



for F{£,L). Hereby we have to take into account that 
not every function F(£, L) that satisfies this equation is a 
physically acceptable solution: an acceptable solution has 
to be non-negative over the entire phase space. In general, 
the problem of solving the integral equation (0) is a de- 
generate problem, because there are in general infinitely 
many distribution functions possible for a given potential- 
density pair. 

Particularly interesting are models for which the dis- 
tribution function and its moments can be computed an- 
alytically. Such models have many useful applications, 
which can roughly be divided into two classes. On the 
one hand, they can improve our general understanding of 
physical processes in galaxies in an elegant way. For ex- 
ample, they can serve as simple galaxy models, in which 
it is easy to generate the starting conditions for A^-body 
or Monte Carlo simulations, or to test new data reduction 
or dynamical modelling techniques. A quick look at the 
overwhelming success of simple analytical models, such 
as the Plummer sphere (Plummer 1911; Dejonghe 1987), 
the isochrone sphere (Henon 1959, 1960), the Jaffe model 
(Jaffe 1983) and the Hernquist model (Hernquist 1990), 
provides enough evidence. On the other hand, analytical 
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models are also useful for the detailed dynamical mod- 
elling of galaxies. For example, in modelling techniques 
such as the QP technique (Dejonghe 1989), a dynamical 
model for an observed galaxy is built up as a linear com- 
bination of components, for each of which the distribution 
function and its moments are known analytically. As a re- 
sult, the distribution function and the moments of the final 
model are also analytical, which obviously has a number 
of advantages. 

Unfortunately, the number of dynamical models for 
which the distribution function is known analytically is 
rather modest. Moreover, most of them consist of dis- 
tribution functions that are isotropic or of the Osipkov- 
Merritt type, and therefore basically one-dimensional. An 
exception is the completely analytical family of anisotropic 
models described by Dejonghe (1987). These models self- 
consistently generate the Plummer potential-density pair, 
a simple yet useful model for systems with a constant den- 
sity core. 

During the last decade, however, it has become clear 
that, at small radii, elliptical galaxies usually have cen- 
tral density profiles that behave as with 0^7^ 2.5 
(Lauer et al. 1995; Gebhardt et al. 1996). Such galaxies 
can obviously not be adequately modelled with a con- 
stant density core. This has stimulated the quest for sim- 
ple potential-density pairs, and corresponding distribution 
functions, with a central density cusp. The first effort to 
construct such models was undertaken by Ciotti (1991) 
and Ciotti & Lanzoni (1997), who discussed the the dy- 
namical structure of stellar systems following the i?^/™ 
law (Sersic 1968), a natural generalization of the empiri- 
cal law of de Vaucouleurs (1948). A major drawback 
of this family, however, is that the spatial density and the 
distribution function can not be written in terms of ele- 
mentary functions (see Mazure & Capelato 2002). A more 
useful family is formed by the so-called 7-models (Dehnen 
1993; Tremaine et al. 1994), characterized by a density 
proportional to r""* at large radii and a divergence in the 
center as with ^ 7 ^ 3. The dynamical structure 
of models with this potential-density pair has been exten- 
sively investigated (e.g. CaroUo, de Zeeuw & van der Marel 
1995; Ciotti 1996; Meza & Zamorano 1997), but only for 
isotropic or Osipkov-Merritt type distribution functions. 
Simple analytical models with a more general anisotropy 
structure are still lacking. 

In this paper we construct a number of families of 
completely analytical anisotropic dynamical models that 
self-consistently generate the Hernquist (1990) potential- 
density pair. It is a special case of the family of 7- 
models, corresponding to 7 = 1. In dimensionless units, 
the Hernquist potential-density pair is given by 

V'(r-) = (2a) 



1 



1 + r 
1 



27r r(l -I- rf 



(2b) 



More precisely, the surface brightness profile has the form 
1 (2 + i?2)X(i?)-3 



^^^^ 27r (l-i?2)2 
with X{R) a continuous function defined as 



(3) 



X{R) 



(1 - i?2)-i/2 arcsechi? 
(i?2 _ l)-i/2 arcseci? 



for s; i? < 1, 
for 1 ^ i? ^ cx). 

(4) 



The paper is organized as follows. The general the- 
ory on the inversion the fundamental integral equation 
(|l]) is resumed in Section 2. Each of the subsequent sec- 
tions is devoted to special cases of this inversion tech- 
nique and the corresponding family of Hernquist mod- 
els. Isotropic models are the most simple ones; Hernquist 
(1990) showed that, in the special case of isotropy, the 
distribution function and its moments can be calculated 
analytically. We repeat the most important characteristics 
of the isotropic Hernquist model in Section 3. In Section 
4 we construct a one-parameter family of models with a 
constant anisotropy. In Section 5, a two-parameter fam- 
ily of Hernquist models is constructed by means of the 
Cuddeford (1991) inversion technique. These models have 
an arbitrary anisotropy in the center and are radially 
anisotropic at large radii. On the contrary, in Section 6, a 
two-parameter family is constructed that has a decreasing 
anisotropy profile, with arbitrary values for the anisotropy 
in the center and the outer halo. Finally, Section 7 sums 
up. 

2. The construction of anisotropic models 

A general discussion on the inversion of the fundamental 
equation ([!]), and hence on the construction anisotropic 
distribution functions for a given spherical potential- 
density pair, is presented by Dejonghe (1986). The key 
ingredient of the inversion procedure is the concept of the 
augmented mass density p(jp,r), which is a function of 
potential and radius, such that the condition 



p{ip{r),r) = p{r) 



(5) 



As the density diverges as l/r for r ^ 0, the surface 
brightness I{R) will diverge logarithmically for i? ^ 0. 



is satisfied. The augmented mass density is in fact equiva- 
lent to the distribution function F{£^ L): with every aug- 
mented density p{tp, we can associate a distribution 
function F{£, L) and vice versa. There exist various tran- 
sition formulae between these two equivalent forms of a 
dynamical model, amongst others a formalism that uses 
combined Laplace-Mellin transforms. 

Besides providing a nice way to generate a distribution 
function for a given potential-density pair, the augmented 
density is also very useful to calculate the moments of 
the distribution function. The anisotropic moments are 
defined as 

M2„.2m(r) =2^// ^(^. L) 1,2™+! dVr dvt, (6) 
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where Vt = y'^'g + is the transverse velocity. One can 

derive a relation that links the higher-order moments to 
the augmented mass density p, when written explicitly as 
a function of ■0 and r, 



X / i^-r) 



2™+" r(n+ i) 
^/tt r(m + n) 
d" 



/\m+n — l 







[dr 



.2\m 



[r2"p(V.',r)]dV''. (7) 



In particular, the radial and transverse velocity disper- 
sions can be found from the density through the relations, 



p{r) 



p(V'',r)d^', 



p{r) 



(8a) 
(8b) 



By means of these functions, we can define the anisotropy 
f3{r) as 



Pir) = 1 



2a2(r) 



(9) 



We will in this paper consider augmented densities which 
are separable functions of -0 and r, and we introduce the 
notation 



p{i^,r) = f(ip)g{r). 



(10) 



For such models, the anisotropy can be directly calculated 
from the augmented density as 



as a result of the formulae (I 



(11) 



), m and (10). 



3. Isotropic models 

3.1. Background 

The simplest dynamical models are those where the aug- 
mented density is a function of the potential only, p — 
piip)- For such models, the distribution function is only a 
function of the binding energy, i.e. the distribution func- 
tion is isotropic. In this case, the integral equation (|l|) can 
be inverted to find the well-known Eddington relation 



m = 



dp dip 



2V27r2 d£ Jo d0 V^^' 



(12) 



For such isotropic models, we do not use the general 
anisotropic moments (^ , but define the isotropic moments 
as 



/i2«(r) =4^ / F{£)v 



,,2n+2 



dv. 



(13) 



Similarly as for the anisotropic moments, we can derive a 
relation that allows to calculate the augmented isotropic 



moments from the augmented density p{ip). Indeed, they 
satisfy the relation (Dejonghe 1986) 



A2«('0) 



(2n+ 1)!! 
(n- 1)!! 



(^_^')"-ip(^')dVy. (14) 



In particular, we obtain an expression for the velocity dis- 
persion profile by setting n = 1, 



<y\r) 



1 



p{r) 



p(V')dV''- 



(15) 



3.2. The isotropic Hernquist model 



The isotropic model that corresponds to the potential- 
density pair (pa|b) is described in full detail by Hernquist 
(1990). We restrict ourselves by resuming the most impor- 
tant results, for a comparison with the anisotropic models 
discussed later in this paper. The augmented density reads 



1 



(16) 



27r 1 - ^ 

Substituting this density into Eddington's formula ( [T2| 
yields the distribution function 



F{£) 



1 



8V27r3 
N/g (1 - 2g) (8g2 - 8g 



3) Sarcsin 

(l-£)5/2 



(17) 



Combining the density (|16| ) with the general formula ( 14 ) 
gives us the moments of the distribution function, 

3-2"+9/2 T{n+\) 



(27r)3/2 T{n + 5) 



V'"+SFi(5,l;n + 5;V). 



(18) 

For all n ^ 0, this expression can be written in terms of 
rational functions and logarithms. For example, for the 
velocity dispersions, we obtain after substitution of the 
Hernquist potential (pa|). 



CT^(r) = r (1 + rf In 



1 



r (25 + 52r -I- 42r2 -t- 12r3) 



12(1 



(19) 

in agreement with equation (10) of Hernquist (1990). From 
an observational point of view, it is very useful to obtain 
an explicit expression for the line-of-sight velocity disper- 
sion. For isotropic models, the line-of-sight dispersion is 
easily found by projecting the second-order moment on 
the plane of the sky, i.e. 



I{R) 



p{r) (j'^{r) r dr 



i?2 



(20) 



For the Hernquist model this yields after some algebra 



I{R)<jI{R) 



1 



247r(l-i?2)3 
3i?2 (20 - 35i?2 + 287?"* 



(6 - 65i?2 + 68i?'* - 24i?6 



R 
"2 ' 



(21) 
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4. Models with a constant anisotropy 

4.1. Background 

A special family of distribution functions that can easily 
be generated using the technique outlined in Section 2 
corresponds to models with a density that depends on r 
only through a factor r~^^, i.e. 

It is well known that such densities correspond to models 
with a constant anisotropy (i.e. Binney & Tremaine 1987), 
which can easily be checked by introducing g{r) = 7--2/3 



in the formula (11). For a given potential-density pair, a 
family of models with constant anisotropy can hence be 
constructed by inverting the potential as r(ip), and defin- 



ing 



/(V') = p(^)(r(^))'''. (23) 

Notice that such models are not the only constant 
anisotropy models corresponding to a given potential- 
density pair, as argued section 1.5 of Dejonghe (1986). 
This family, however, is very attractive due to its relative 
simpleness. In particular, the corresponding distribution 
function is a power law in L, and can be found through 
an Eddington-like formula, 

^(^' ^) ^ r(l-/3)r(i+/3) 

" d/ dip 



d£ Jo dip {£-ip)U2'i3- (2"^) 

For the moments of the distribution function, the relation 
(^ can be simplified to 

2"'+" r(n+l)r(m+l~/?) 

M2«,2m W-, r) = — = — ■ — — — -— r 

Vtt r{m + n)r(l - P) 

X r(V-^')™+"~'/(^')dV''- (25) 



(26) 



In particular, the radial velocity dispersions reads 
1 /■■^('') 



fmr)) Jo 



/(V')d^'. 



4.2. Hernquist models with a constant anisotropy 
4.2.1. The distribution function 

Applying the formula ( p3|) to the Hernquist potential- 
density pair (pa|b) yields 

f{i^) = ^^'-''{^-^f'-'- (27) 

Substituting this expression into the general formula ( p4[ ) 
gives us the corresponding distribution function 



F{£,L) 



2" 



r(5 - 2/3) 



(27r)5/2 Y(l-p)Ta-p) 



gi-P U - 2/3, 1 - 2/3; ^ - /3; f ) . (28) 



This expression reduces to the isotropic distribution func- 
tion (|^ for /3 = 0, as required. Whether the expression 
( p8| ) corresponds to a physically acceptable distribution 
function for a given value of /3 depends on the condition 
that the distribution function has to be positive over the 
entire phase space. 

It is no surprise that the distribution function is not 
positive for the highest possible values of because mod- 
els where only the radial orbits are populated can only 
be supported by a density profile that diverges as r~'^ 
or steeper in the center (Richstone & Tremaine 1984). It 
turns out that the distribution function (|2^) is everywhere 
non-negative for /3 ^ ^. 

For all integer and half-integer values of /3, the hy- 
pergeometric series in ( p8| ) can be expressed in terms of 
elementary functions. Very useful are half-integer values of 
/3, because the energy-dependent part of the distribution 
function can then be written as a rational function of £. 
For integer values of /3, the hypergeometric series can be 
written as a function containing integer and half-integer 
powers of £ and \ — £ and a factor arcsin similar to 
the isotropic distribution function (p7|). 

The limiting model [3 — ^is particularly simple. It has 
an augmented density that is a power law of potential and 
radius. 



p(V',r) 



(29) 



and the corresponding distribution function simply reads 
on 1973), 



O Cl 



(30) 



This model is a special case of generalized polytropes dis- 
cussed by Fricke (1951) and Henon (1973). 

In Fig. we compare the distribution functions of the 
radial model with /5 = ^ and the tangential model with 
/3 = — 2 with the distribution function of the isotropic 
Hernquist model. The distribution functions are shown 
by means of their isoprobability contours in turning point 
space, which can easily be interpreted in terms of orbits. 
Compared to the isotropic model, the radial model prefers 
orbits on the upper left side of the diagram, with an apoc- 
enter much larger than the pericenter, i.e. elongated or- 
bits. The isoprobability contours of tangential models on 
the other hand lean towards the diagonal axes where peri- 
center and apocenter are equal, i.e. nearly-circular orbits 
are preferred. 

4.2.2. The velocity dispersions 

By means of substituting the expression (|2^) into the gen- 
eral formula (Q), we can derive an analytical expression for 
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Fig. 1. The distribution function of the Hernquist models with a constant anisotropy, represented as isoprobability 
contours in turning point space. The distribution functions in solid lines represent a radial model with (3 — ^ (left 
panel) and a tangential model with f3 — —2 (right panel). The dotted contour lines in both panels correspond to the 
isotropic Hernquist model. 



0.6 



0.4 - 



0.2 - 



- 
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Fig. 2. The velocity dispersion of the Hernquist models 
with a constant anisotropy. The upper and lower panels 
show the radial velocity dispersions ar{r) and the line-of- 
sight velocity dispersion ap{R). The profiles are shown for 
different values of the anisotropy parameter /3: plotted are 



f3 = 



1 1 

2 ' 4 ' 2 ' 



0, 



-2, —5 and the limit case (3 —> — oo. 



We are mainly interested in the velocity dispersions, which 
can be conveniently written by means of the incomplete 
Beta function (Abramowitz & Stegun 1972), 

aUr) = (i + rf (s - 2/3, 2/?) . (32) 

For all anisotropics (3 < ^, the radial dispersion equals 
zero in the center of the galaxy, rises until a maximum 
and then decreases again towards zero for r —^ oo. The 
asymptotic behavior for r 3> 1 is 



all moments of the distribution function, 

^ ^ ^n+m-i r(5 - 2/3) r(n + I) r(m + 1 - (3) 
M2„,2m(^,r) - ^^75 r(m + n + 5- 2/3) r(l-/3) 

X 2i^i (5-2/3, l-2/3;m + n + 5-2/3; V). (31) 



1 1 



5-2/3 r 



(33) 



The expression (|32| ) can be written in terms of elementary 
functions for all /3 with 4/3 an integer. For the integer and 
half-integer values of /3, the expression involves polynomi- 
als in r and a factor ln(l -I- 1/r), very analogous to the 
expression ( p^ of the isotropic Hernquist model. For the 
quarter-integer values of /3, it contains polynomials and 
square roots in r and a factor arccotg ^/r. 

Particular cases are the models that correspond to the 
most radial and tangential distribution functions. On the 
one hand, the limit case (3 = ^ has the simple velocity 
dispersion profiles 

1 1 
4 1 + r' 



:(r) 



(34) 



Particular about this dispersion profile is that it assumes a 
finite value in the center. On the other side of the range for 
possible anisotropies, we can consider the limit /3 — > — oo, 
which corresponds to a model with purely circular orbits. 
For such a model, the radial dispersion is of course identi- 
cally zero, whereas the transverse velocity dispersion is the 
circular velocity corresponding to the Hernquist potential. 



1 



2 (1 



(35) 
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In the top panel of Fig. || we plot the radial veloc- 
ity dispersion profile for various models with a different 
anisotropy (3. Both at small and large radii, the radial 
dispersion is a decreasing function of /3, as expected. 

The line-of-sight velocity dispersion for anisotropic 
models is found through the formula 



P{r) (J^s^r, R)rdr 



(36) 



where ctios (T^ R) is the velocity dispersion at the position 
r on the line of sight R in the direction of observer. It is 
a linear combination of the radial and transverse velocity 
dispersions in this point, 



aUr,R)= 1- 



R^ 



al{r) 



R^ 

2r2 



We can equivalently write 



aliR) 



liR) 



n2 



p{r) cr^(r) r dr 



i?2 



(37) 



(38) 



For the general Hernquist models with a constant 
anisotropy, the integration (|3^) cannot be performed ana- 
lytically. But for all integer and half-integer /3's, <Jp{R) can 
be expressed in terms of polynomials and the function X, 
defined in equation (^). For example, for the limit model 
/3 = i we obtain after some algebra 



IiR)4{R) 



487r(l-i?2)3 
3(4 - 14i?2 + 35i?4 _ 28R^ + 8R^) X{R) 

- (28 - 57i?2 + eSi?" - 24i?6) +-. (39) 



For the other limit model, the one with only circular or- 
bits, we find 



I{R)al{R) 



i?2 



' 487r(l-i?2)4 
-(120 - 120i?2 + 189i?'* - lOSi?^ + 24i?^) X{R) 



(154 - 117i?^ + %2R^ - 2AR^) 



R 
I' 



(40) 



in agreement with Hernquist (1990). 

In the bottom panel of Fig. ^ we plot the line-of-sight 
dispersion profiles for a number of different values of (3. 
The behavior of the individual profiles is analogous to the 
spatial dispersion profiles: except for the /3 = ^ model, 
which has a finite central dispersion, the profiles start 
at zero in the center, rise strongly until a certain maxi- 
mum and then decrease smoothly towards zero at large 
projected radii. The behavior for i? ^ 1 can be quanti- 
fied if we introduce the asymptotic expansion ( |33| ) into the 
formula iM, 



4/3 



IStt V5-2/3 



1 

R 



(41) 



The dependence of the line-of-sight dispersion as a func- 
tion of the anisotropy is depends strongly on the projected 
radius: at small projected radii, (Jp decreases with increas- 
ing /3, whereas for large radii, dp increases with increas- 
ing (3. In other words, radial models have a larger central 
and a smaller outer line-of-sight dispersion than their tan- 
gential counterparts. This is a direct consequence of the 
weight of the radial and transversal velocity dispersions in 
the linear combination (|37|): at small projected radii, the 
radial dispersion contributes the dominant term, whereas 
for the outer lines of sight, the transversal dispersion term 
dominates. 

5. Models with increasing anisotropy 

5.1. Background 

Osipkov (1979) and Merritt (1985) developed an inversion 
technique for a special class of distribution functions that 
only depend on energy and angular momentum through 
the combination 



2rV 



(42) 



with Ta the so-called anisotropy radius, with the additional 
condition that 



F{£,L) = 



for g < 0. 



(43) 



Such models correspond to an augmented density of the 
form 



(44) 



p(V,r)= 1 + ^ f{^). 



In this case, the fundamental integral equation (|i|) can be 
inverted in a similar way as the Eddington relation. 



F{£,L) 



1 



d/ dV 



2\/27r2 dg 7o dV- VQ^' 



(45) 



The anisotropy /3(r) for the Osipkov-Merritt models can 
be found by means of formula (fl]). 



r.2 • 



(46) 



These models are hence isotropic in the center and com- 
pletely radially anisotropic in the outer regions. The pa- 
rameter ra determines how soon the anisotropy turns from 
isotropic to radial. In particular, for r^, — > oo, g is noth- 
ing else than the binding energy, and the Osipkov-Merritt 
models reduce to the isotropic models. 

The Osipkov-Merritt models were generalized by 
Cuddeford (1991), who considered models which corre- 
spond to an augmented density of the form 



p(-0,r) 



-2/3o 



1 - 

r'i 



(47) 
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These models reduce to the Osipkov- Merritt models if we 
set /3o = 0. Within this formalism, the distribution func- 
tion can be calculated in a similar way as for the Osipkov- 
Merritt models. The distribution functions have the gen- 
eral form 

Fi£,L) = Fo{Q)L-^'''>, (48) 

with the additional condition (|4^). The solution of the 
integral equation (|l|), valid for f3o < 1, reads in this case 



Fi£,L) 



where 



2/3o 



(27r)3/2 r(l-/3o)r(l-0) 



dg 



^3 d"/ d^A 



dV'™ (Q - ipy 



m = 1 -f- int ( - - /3o 



= frac Po 



r, (49a) 

(49b) 
(49c) 



The most interesting cases are those where /3o is either 
integer or half-integer. For integer values of /3o, the general 
formula (49a be) reduces to 



F{£,L) 



2^2 7r2 r(l - /?n) dQ J„ d^i-Z^o 



(50) 



For half- integer values of /3o, the integral equation (|l|) is a 
degenerate integral equation, which can be solved without 
a single integration (Dejonghe 1986; Cuddeford 1991)]^ 



F{£,L)^ 



2/3o 



(2^)3/2 r(l - f3o) 



d^~P"/ 
dipi 



-00 



(51) 



As for the Osipkov-Merritt models, the anisotropy of the 
Cuddeford models has a very simple functional form, 
which can be found through (pd]). 



(52) 



They hence have an anisotropy Po in the center]^ and be- 
come completely radially anisotropic in the outer regions. 
The anisotropy radius is again a degree for how quick 
this transition takes place. In particular, for ^ cxd, 
the Cuddeford models reduce to models with a constant 
anisotropy P — Pq. Because the range of values for Po for 



which the inversion (49a be) is mathematically defined is 
only restricted by /3o < 1, distribution functions can in 
principle be calculated with any degree of anisotropy in 
the center, ranging from very radial to extremely tangen- 
tial. Whether these distribution functions correspond to 
physically acceptable solutions depends on the positivity, 
however. 



^ Apparently, a factor r( | — /3o) is missing in the denominator 
in formula (30) of Cuddeford (1991). 

^ This explains why we prefer the parameter f3o above a = 
— /3o originally adopted by Cuddeford (1991). 



5.2. Hernquist models with increasing anisotropy 
5.2.1. The distribution function 



For the Hernquist potential-density pair (|2a|b), the aug- 
mented density corresponding to the Cuddeford formalism 
is readily calculated. We obtain 

n-l 1-00 



1 + A 



1 - V 



i-200 



(1 - iijY-'^f^« ' 



(53) 



where we have set A = l/r^. Combining this expression 



with the general Cuddeford solution (49a be) we can ob- 
tain distribution functions that self-consistently generate 
the Hernquist potential-density pair, and which have an 
arbitrary anisotropy in the center and a completely radial 
structure in the outer regions. In order to represent phys- 
ically acceptable dynamical models, it is necessary that 
these distribution functions are positive over the entire 
phase space, i.e. F{£,L) ^ for ^ Q ^ 1. Before try- 
ing to actually calculate the distribution functions, it is 
useful to investigate which region in the (/3o,A) parame- 
ter space corresponds to physically acceptable distribution 
functions. 

First of all, it is obvious that the models with Pq > ^ 
will not correspond to non-negative distribution functions: 
the distribution function is already too radial for A = 
(Sec. 4.2.1), and will become even more radial for larger 
A. We can therefore limit the subsequent discussion to 
Po ^ ^. Now consider such a fixed value Pq, and con- 
sider all Cuddeford models corresponding to this central 
anisotropy. For A — *■ 0, the Cuddeford model reduces to 
the model with constant anisotropy Pq, which is physi- 
cally acceptable (Sec. 4.2.1). For A oo, the distribution 
function will only consist of radial orbits, for which the 
distribution function is not positive. It can therefore be 
expected that, for a given value of /3o ^ ^, a range of A's 
is allowed, starting from up to a certain Amax- 

Next, we have to investigate how A,nax varies with Pq, 
i.e. which anisotropy radii are allowed for a given central 
anisotropy ? Distribution functions with a strong central 
tangential anisotropy and a small anisotropy radius are 
likely to be negative. Indeed, consider the orbital struc- 
ture of such a galaxy. Because the outer regions of the 
galaxy (r ^ Tq) are strongly radially anisotropic, the vast 
majority of the stars there must be on nearly radial orbits. 
These stars also pass through the central regions, where 
they will contribute to the central density and radial ve- 
locity dispersion as well. The smaller the value of ra, i.e. 
the larger the value of A, the stronger the contribution 
of stars on such nearly radial orbits. In order to create 
a core where the anisotropy is tangential, a large num- 
ber of stars hence have to be added which move on tightly 
bound nearly circular orbits. But we are limited from keep- 
ing on adding such stars, because we cannot exceed the 
spatial density of the Hernquist profile, which has only a 
fairly weak divergence. We therefore expect that no 
Cuddeford models will exist beyond a certain minimal Pq 
(except for the degenerate case of the constant anisotropy 
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Table 1. The range of anisotropy radii which give rise 
to a positive distribution function of the Cuddeford type, 
consistent with the Hernquist potential-density pair. For a 
given value of Pq, this range corresponds to ^ A ^ A„iax, 
or equivalently, to Ta.min ^ ^ oo. 









^ -1.500 


0.000 


OO 


-1.375 


1.764 


0.753 


-1.250 


3.598 


0.527 


-1.125 


5.550 


0.424 


-1.000 


7.582 


0.363 


-0.875 


9.680 


0.321 


—0.750 


11.83 


0.291 


-0.625 


14.02 


0.267 


-0.500 


16.23 


0.248 


-0.375 


18.51 


0.232 


-0.250 


20.57 


0.220 


-0.125 


22.61 


0.210 


0.000 


24.42 


0.202 


0.125 


25.87 


0.197 


0.250 


26.70 


0.194 


0.375 


26.42 


0.195 


0.500 


24.00 


0.204 




Fig. 3. The region in (/3o, A) space corresponding to a pos- 
itive distribution function of the Cuddeford type, consis- 
tent with the Hernquist potential-density pair. 

models, which have no radial anisotropy at large radii). 
Moreover, it can be expected that for models with a tan- 
gential central anisotropy, the range of anisotropy radii is 
more restricted than for models with a radial or isotropic 
central anisotropy, i.e. that Amax(/?o) is a increasing func- 
tion of i3q. 

By numerical evaluation of the integral in equation 
(49a), we calculated Amax(/3o) for a number of values for 
Po (Table |l|). The region in parameter space where the 
Cuddeford-Hernquist models are physical is shown in Fig. 
^. Notice that all models with /3o ^ — | and A > are neg- 
ative at some point in phase space and are thus unphys- 
ical: the Hernquist potential-density pair can support no 
(non-degenerate) distribution functions of the Cuddeford 
type with a central anisotropy /3o =^ — § • 

We are primarily interested in those models where the 
distribution function can be expressed in terms of elemen- 
tary functions. This is of course possible for all half-integer 
values of (3q, because the calculation of the distribution 



function involves no integrations. Also for the integer val- 
ues of /So, the distribution function can be calculated ana- 
lytically, through the formula ( pO| ) . Because of the limited 
region in (/3o, A) space where Cuddeford models are non- 
negative, this leaves us with four models with analytical 
distribution functions, corresponding to Pq = ^, 0, — ^ 
and —1. The most simple of them is the case /3o = 5, for 
which we obtain 



A(3Q2_5Q + 2) 



1 + A 



(54) 



For this distribution function, it is straightforward to 
check that it remains positive for ^ A ^ 24, in agree- 
ment with the numerical result in Table ^. For Pq — 0, we 
recover the Osipkov-Merritt model, 



Fi£,L) 



8%/27r3 

+Vo(i- 



3 arcsin ^/Q 
2Q) 



Q)5/2 



8A 



(55) 



in agreement with Hernquist (1990). For the two other 
cases, Po — and Po = —1, the distribution function 
can also be written in terms of elementary functions, but 
the expressions are somewhat more elaborate. 

In Fig. I we show the distribution function of the 
Cuddeford type for four different models. The models on 
the top row have a radial central anisotropy, whereas those 
in the bottom panels have a tangential anisotropy in the 
center. The left and right column correspond to two dif- 
ferent values of the anisotropy radius. The dotted distri- 
bution functions on the background are the distribution 
functions with a constant anisotropy Pq. 

The character of the Cuddeford models can directly be 
interpreted from these figures. Compared to the constant 
anisotropy models, the Cuddeford models have a much 
larger fraction of stars on radial orbits, visible for both 
models with radial and tangential central anisotropy. The 
most conspicuous feature of each of the Cuddeford dis- 
tribution functions is that the right part of the (r_,r_(-) 
diagram is completely empty, i.e. at large radii only the 
most radial orbits are populated, which is necessary to 
sustain the radial anisotropy. The boundary of the region 
in turning point space beyond which no orbits are popu- 
lated can be calculated by translating the equation Q = 
in terms of the turning points r_ and r+. 



1 + Xr'i 1 + A ' 



(56) 



When we substitute the Hernquist potential (pa|), we can 
actually calculate the range of allowed orbits, 

^ r_ < ?-c,max (57) 

(1 + r_) + Jl + 2r_ +r'i+ 4Ari 
r- r+ L_ , (58) 
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Fig. 4. Comparison of the distribution function corresponding to Hernquist models of the Cuddeford type and 
Hernquist models with a constant anisotropy. The distribution functions are represented as isoprobability contours 
in turning point space. The solid lines correspond to the Cuddeford distribution functions, with the parameters f3o 
and A displayed in the down left corner of each diagram. The dotted lines represent the distribution functions of the 
corresponding Hernquist models with a constant anisotropy /No- 



where Tcmax represents the radius of the largest allowed 
circular orbit for a given A, 

32/3 VA (971+ V8TA^)V3- ^''^ 

Obviously, the larger A, the more restricted the range of 
allowed orbits, because the transition to radial anisotropy 
occurs at smaller radii for large values of A. This can be 
seen when comparing the left and right panels of Fig. 0. 



5.2.2. The velocity dispersions 

In order to calculate the radial velocity dispersion associ- 
ated with models of the Cuddeford type, we use the gen- 
eral formula (pa|). After some manipulation, we obtain 



(1 + Ar2)i-Ai 



{1 + Xr'^)^~l^»dr' 
j,/l-2/3o(l -I- r')^ 



(60) 



In general, this integral needs to be evaluated numerically, 
but for the four models with integer and half-integer values 
of /3o, it can be performed analytically. For example, for 
the Osipkov-Merritt model Pq — 0, we find 



r (1 + r)" 



In 



1 



1 + Ar^ \ r 
r (25 -I- 52r -I- 42r'^ + 12r^) - A r (1 + 4r) 
12(l + r)(l + Ar2) ' 



(61) 



which reduces to the isotropic dispersion (19) for A — 0. 
For the other integer and half-integer values of /3q, the ra- 
dial dispersion can also be expressed in terms of algebraic 
functions and logarithms, but the expressions are some- 
what more elaborate. 

In the top panels of Fig. || we plot the radial veloc- 
ity dispersion profiles for Hernquist-Cuddeford models, for 
varying /3o and varying A (left and right panels respec- 
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Fig. 5. The radial (upper panels) and line-of-sight (lower panels) velocity dispersion profiles of the Hernquist-Cuddeford 
models. The different curves in the two left panels correspond to models with the same anisotropy radius A = 1 (i.e. 
Ta = 1), but with a different central anisotropy parameter Pq: plotted are /3 = ^, 0, — ^ and —1. The dotted curves 
are the dispersion profiles of the (hypothetical) completely radial Hernquist model, which corresponds to /3o = 1- The 
two panels on the right-hand side contain the dispersion profiles of models with the same central, slightly tangential, 
anisotropy parameter /3o = —5, but with a varying anisotropy radius. The various curves correspond to A = (black 
solid line), 0.2, 1, 3, 8 and A,nax = 16.23. Again, the dotted curves are the dispersion profiles of the (hypothetical) 
completely radial Hernquist model, which corresponds to A — > 00. 



lively). The behavior of function of /3o is pre- 

dictable. At small radii, the different models have a dif- 
ferent behavior, with the largest dispersion for the most 
centrally radial models. At large radii they all have a sim- 
ilar, purely radial, orbital structure, and as a consequence 
their dispersion profiles all converge towards a single pro- 
file. This limiting profile is the radial velocity dispersion 
profile that corresponds to the (hypothetical) model with 
a completely radial orbital structure, which we can obtain 
by either setting /3 = 1 in the expression (|3^), or setting 
/3o = 1 in the expression (60), 



1 1 + 4r 



12 r(l + r) 



(62) 



For a fixed central anisotropy, the behavior of the radial 
dispersion as a function of the anisotropy radius also fol- 
lows a simple trend: the Or profiles increase with increasing 
A, and the curves are all bounded by two limiting profiles: 
on the one hand the dispersion profile ( p2[ ) of the con- 
stant anisotropy model (obtained by setting A = 0), and 
on the other hand the hypothetical dispersion profile ( |62[ ) 
of the purely radial model (which corresponds to A ^ cxd). 
Dispersion profiles with large A will more quickly lean to- 
wards the purely radial profile than models with small 
A, because the transition to a strongly radial anisotropy 
occurs zi r ^ 1/VA. 



The bottom panels of Fig. || show the linc-of-sight 
velocity dispersion of the Hernquist-Cuddeford models. 
These profiles had to be calculated numerically. The de- 
pendence of the line-of-sight dispersion upon A and /3o can 
be easily interpreted. In particular, the line-of-sight dis- 
persion profiles of the Cuddeford models tend towards the 
line-of-sight dispersion profile of the hypothetical purely 
radial Hernquist model, which reads 



R 



I{R)cT;iR) = - 



1 



1 



^ 96i? 487r (1 - i?2)2 
i?2 (20 - 29i?2 + 12^^^) X{R) + (2 - 7R^ 



6. Models with decreasing anisotropy 



(63) 



6.1. Background 

In order to construct dynamical models with a decreas- 
ing anisotropy, i.e. with a tangentially anisotropic halo, 
no special inversion techniques exist, such that we have to 
rely on the general formulae of Dejonghe (1986) to invert 
the fundamental integral equation (^. A disadvantage is 
that these formulae are numerically unstable. Their useful- 
ness is therefore actually restricted to analytical models. 
But this is not straightforward: a direct application of the 
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inversion formulae to an arbitrary analytical augmented 
density p{'ip,r), even if its looks rather simple, can result 
in a cumbersome mathematical exercise, because the in- 
version formulae are quite elaborate. 

A useful strategy to construct models with a tangential 
halo without the need to cope with the complicated gen- 
eral formulae, is to profit from the linearity of the integral 
equation (^. In particular, it is very interesting to gener- 
ate augmented densities p('!/',r), which can be expanded 
in a series of functions pk{tp,r), which depend on r only 
through a power law. 



(64) 



Each of the augmented densities pk{ip,r) corresponds 
to a dynamical model with a constant anisotropy (3k- 
Combining the lin eari ty of the integral equation (|l|) with 
the results of Sec. O , we find that the distribution func- 
tion corresponding to the density (joj) reads 

F(£,L) =^^^fc(£,L), (65a) 



with 



2/3. 



(2^)3/2 r(i-/5fc)r(i + /3fc) 



dV- 



d£ 



dh 

dV- {£ -ipy/^'P"' 



(65b) 



Equivalently, the moments of the distribution function can 
be derived from the series expansion. 

6.2. Hernquist models with decreasing anisotropy 

6.2.1. Construction of the density function 

For every potential ?A(r), we can create an infinite number 
of functions Z{il;, r) which satisfy the identity Z{ip{r), r) = 
1. For the Hernquist potential, we can easily create such 
a one-parameter family of functions Zn^^p, r), 



1, 



(66) 



with n a natural number. If we multiply this family 
with the density function ( p7| ) of the constant anisotropy 
Hernquist models, we create a new two-parameter family 
of dynamical models, that will self-consistently generate 
the Hernquist potential-density pair, 



1 



4-2Al+n 



(1 

2tt (1 - tpy-^Po r- 



200 



(67) 



Defining a new parameter /3oo = /3o ~ we can write this 
augmented density also as 



p{ij,r) 



1 



4-2/3 o, 



(1 



n2(/3o-/3oo) 



2tt (1 - tpy-^i^o 



r2Po 



(68) 
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Fig. 6. Comparison of the distribution function corre- 
sponding to Hernquist model with increasing anisotropy 
with f3o = ^ and f3oo = and the constant anisotropy 
model with f3 = —2. The distribution functions are repre- 
sented as isoprobability contours in turning point space. 

Because we assumed that n is a natural number, we can 
expand the binomial in the nominator of the density ( |67[ ) , 
and write it in the form (p4), with 



2^ \k) (1-^)1-2/30 



(69a) 
(69b) 



with ^ k ^ n. The reason why we chose Pq and Poo as 
parameters becomes clear when we look at the expression 
for the anisotropy corresponding to this family of den- 
sity functions - for the moment being without bothering 
whether the density corresponds to a physically accept- 
able distribution function. By means of the formula (pi]), 
we obtain 



Pir) 



pQ+ Poor 
1 + r 



(70) 



The anisotropy equals Po in the center and increases to 
Pea at large radii. Because n can in principle assume any 
natural number, this family of density functions hence cor- 
responds to dynamical models which can grow arbitrarily 
tangential in the outer regions. In particular, by setting 
n = we recover the models with constant anisotropy 
f3^/3o=Poo from Sec. O. 



6.2.2. The distribution function 

We can calculate the dist ribu tion function of these models 
by applying the recipe ( |65b| ) to each of the components 
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9a|b). We obtain after some algebra 

2/3o 



F{£,L) = 



E 



(27r)5/; 



1 



7_k 
2^2 



fe=0 

X 2F1 ( 5-2/3oo,l-2/3o 



2/3oo+/3o;fV (71) 



This family of models is restricted by the condition (3q ^ 
^ , because for higher values of /Jq the distribution function 
becomes negative. For all half-integer and integer values 
of Pq (and therefore also of (3oo), this expression can be 
written in terms of elementary functions, very analogous 
with the distribution functions of the constant anisotropy 
models: the expression contains integer and half-integer 
powers of £ and 1 — £ and a factor arcsin y/£. The models 
characterized by /3o = ^ are of a particular kind. The 
hypergeometric functions in equation ( [7l| ) disappears for 
/3o = 5 , such that the distribution function can be written 
as a finite power series of V£ and L. 

An interesting characteristic of these models is re- 
vealed when we look at the asymptotic behavior of the 
distribution function at large radii, i.e. for £ 0. The 
term corresponding to k — n will contribute the dominant 
term in the sum ([Tw, such that we obtain 



Fi£,L) 



2/3. 



r(5 - 2/3, 



(27r)5/2 r(l-/3^)r(|-/3oo) 

X £-5/2-/3^ 



(72) 



This expansion is at first order independent of f3o, such 
that all models with the same /3oo will have a similar be- 
havior at large radii. In particular, all distribution func- 
tions corresponding to a particular /3oo will at large radii 
behave as the Hernquist model with constant anisotropy 
(3 = Poo- This is illustrated in Fig. ^, where we compare 
the distribution function of a model with a radial core 
and a tangential halo with the constant anisotropy model 
that has the corresponding tangential anisotropy. At small 
radii, the difference between both distribution functions is 
obvious: the former one has more stars on radial orbits, 
whereas the latter prefers to populate circular-like orbits. 
At large radii, however, the isoprobability contours of both 
models agree very well. 

Finally, notice that there is no analogue for this be- 
havior at small radii: not all models with a fixed (3q will 
have a similar behavior for £ — > 1, i.e. at small radii. 

6.2.3. The velocity dispersions 

In order to calculate the velocity dispersion profiles of the 
models of this type, we have various possibilities. W e can 
either calculate the dispersion for each of the n terms ( |69a| ) 
through formula ( p6| ) and sum the results, or directly ap- 
ply the general recipe (|8^b) on the expression (E^). In 
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Fig. 7. The radial (upper panel) and line-of-sight (lower 
panel) velocity dispersion profiles of Hernquist models 
with a decreasing anisotropy. All models have the same 
tangential outer anisotropy Poc = —2, but they have 
a different central anisotropy parameter Pq: plotted are 
Po = 5, 0, — i, —1 and —2 (black line). 



either case, we obtain an expression very akin to the ex- 
pression ( |3^ ) of the models with constant anisotropy, 

aUr) = ri-2A, (1 + ^)3+2^^-2A, (5 - 2/3^, 2/3o) • 

(73) 

This expression can be written in terms of elementary 
functions for all Po with 4/3o an integer (and hence also 
4/3oo an integer). 

Not as a surprise, the asymptotic expressions for (T^{r) 
for r 1 read 

"'"•>»^^-^^ 

i.e. they are similar to the corresponding expansions of the 
constant anisotropy models with P — Poo . This behavior is 
illustrated in the upper panel of Fig. where we plot the 
radial velocity dispersion profile for a set of models with 
varying Pq and a fixed Poo- At small radii, the models have 
different profiles (those with the most radial anisotropy 
have the largest values of ar), but at large radii, they all 
converge towards a common asymptotic expansion. 

The calculation of the line-of-sight velocity dispersion 
is also similar to the case of constant anisotropy. It is 
found that ap{R) can be written in terms of elementary 
functions for all integer and half- integer values of Po, and 
that the asymptotic behavior for i? ^ 1 reads 



2Po 



1 
R 
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which is at first order independent of /3o. An illustration 
is shown in the bottom panel of Fig. |7[ 

7. Conclusions 

Three new families of anisotropic dynamical models 
have been presented that self-consistently generate the 
Hernquist potential-density pair. For all models, in par- 
ticular for the Cuddeford models of Section we checked 
the conditions on the adopted parameters such that the 
distribution is positive, and hence physically acceptable, 
in phase space. 

They host a wide variety of orbital structures: in gen- 
eral, the models presented can have an arbitrary central 
anisotropy, and a outer halo with the same anisotropy, a 
purely radial orbital structure, or an arbitrary, but more 
tangential, anisotropy. In order to produce models that 
have an arbitrary anisotropy in the central regions, and 
a more radial, but not purely radial, anisotropy at large 
radii, the most cost-effective way seems to construct a lin- 
ear combination of a number of 'component' dynamical 
models, such as the ones presented here. This technique 
has been adopted for several years in the QP formalism 
(Dejonghe 1989, for an overview see Dejonghe et al. 2001), 
where most of the components in the program libraries 
have an intrinsically tangential orbital structure. 

For all of the presented models, we have analytical ex- 
pressions for the distribution function and the velocity dis- 
persions in terms of elementary functions. They are hence 
ideal tools for a wide range of applications, for example 
to generate the initial conditions for iV-body or Monte 
Carlo simulations. At this point, a number of remarks are 
appropriate. 

First, very few elliptical galaxies are perfectly spher- 
ical; actually, various observational and theoretical evi- 
dence suggests that many elliptical galaxies are at least 
moderately triaxial (Dubinski & Carlberg 1991; Hernquist 
1993; Tremblay & Merritt 1995; Bak & Statler 2000). 
Unfortunately, an extension of the presented techniques 
to construct analytical axisymmetric or triaxial systems is 
not obvious, because the internal dynamics of such stellar 
systems is much more complicated than in the spherical 
case. Nevertheless, our models can be used as a onset to 
construct numerical axisymmetric of triaxial distribution 
functions with different internal dynamical structures, for 
example by the adiabatic squeezing technique presented 
by HoUey-Bockelmann et al. (2001). 

Second, the models presented here are self-consistent 
models, whereas it is nowadays believed that most ellip- 
tical galaxies contain dark matter, either in the form of 
a central black hole (Merritt & Ferrarese 2001 and ref- 
erences therein) and/or a dark halo (Kronawitter et al. 
2000; Magorrian & Ballantyne 2001). When constructing 
dynamical models with dark matter, an extra component 
must be added to the gravitational potential. For example, 
Ciotti (1996) constructed analytical two-component mod- 
els in which both the stellar and dark matter components 
have a Hernquist density profile and an Osipkov-Merritt 



type distribution function. The models presented in this 
paper can also be extended to contain a dark halo or a cen- 
tral black hole. Indeed, the adopted inversion techniques 
are perfectly suitable for this, because the augmented den- 
sity functions p{il), r) do not necessarily need to satisfy the 
self-consistency condition (^). Adding an extra term to 
the potential does not conceptually change the character 
of the inversion, but it might complicate the mathematical 
exercise. 

Third, we have not discussed stability issues for the 
presented models. The study of the stability of anisotropic 
stellar systems is difficult, and a satisfactory criterion can 
not easily be given. For stability against radial perturba- 
tions, we can apply the sufficient criterions of Antonov 
(1962) or Doremus & Feix (1973), but numerical simu- 
lations have shown that these criteria are rather crude 
(Dejonghe & Merritt 1988; Meza & Zamorano 1997). 
Moreover, the only instability that is thought to be effec- 
tive in realistic galaxies is the so-called radial orbit insta- 
bility, an instability that drives galaxies with a large num- 
ber of radial orbits to forming a bar (Henon 1973; Palmer 
& Papaloizou 1987; Cincotta, Nunez & Muzzio 1996). The 
behavior of galaxy models against perturbations of this 
kind can only be tested with detailed A^-body simulations 
or numerical linear stability analysis. Meza & Zamorano 
(1997) used A^-body simulations to investigate the ra- 
dial orbit instability for a number of spherical models of 
the Osipkov-Merritt type, including the Hernquist model. 
They found that the models are unstable for < 1, which 
significantly restricts the set of models that correspond to 
positive distribution functions (see Table |l|) . It would be 
interesting to extend this investigation to the three fami- 
lies of Hernquist models presented in this paper, but this 
falls beyond the scope of this paper. 

Acknowledgements. The authors are grateful to Andres Meza 
for a careful check on the formulae derived in this paper. 
Fortran codes to evaluate the internal and projected dynamics 
of the presented models are available from the authors. 

References 

Abramowitz M., Stegun I. A., 1972, Handbook of 
Mathematical Functions (Dover Publication Inc., New 
York) 

Antonov V. A., 1962, Vestnik Leningrad University, 19, 96 
(trans, in Princeton Plasma Physics Lab. Kept. PPI- Trans- 
1) 

Bak J., Statler T., 2000, AJ, 120, 110 

Binney J., Tremaine S., 1987, Galaxy Dynamics (Princeton 

University Press, Princeton) 
CaroUo C. M., de Zeeuw P. T., van der Marel R. P., 1995, 

MNRAS, 276, 1131 
Ciotti L., 1991, A&A, 249, 99 
Ciotti L., 1996, ApJ, 471, 68 
Ciotti L., Lanzoni B., 1997, A&A, 321, 724 
Cincotta P. M., Nunez J. A., Muzzio J. C, 1996, ApJ, 456, 274 
Cuddeford P., 1991, MNRAS, 253, 414 
Dejonghe H., 1986, Phys. Rep., 133, 217 
Dejonghe H., 1987, MNRAS, 224, 13 



14 



M. Baes & H. Dejonghe: The Hernquist model revisited 



Dejonghe H., 1989, ApJ, 343, 113 

Dejonghe H., Merritt D., 1988, ApJ, 328, 93 

Dejonghe H., Baos M., De Bruync V., De Rijcke S., Van 
Caelenberg K., Van Hese E., 2001, in Stellar Dynamics: 
from Classic to Modern, eds. L. P. Osipkov & 1. 1. Nikiforov 
(Saint-Petersburg State University), p. 292 

Dehnen W., 1993, MNRAS, 265, 250 

de Vaucouleurs G., 1948, Ann. d'Astrophys., 11, 247 

Doremus J. P., Feix M. R., 1973, A&A, 29, 401 

Dubinski J., Carlbcrg R., 1991, ApJ, 378, 496 

Eddington A. S., 1916, MNRAS, 76, 572 

Pricke W., 1951, Astron. Nachr., 280, 193 

Gebhardt K., Richstonc D., Ajhar E. A., Laucr T. R., Byun 
Y.-L, Kormcndy J., Dressier A., Fabcr S. M., Grillmair C., 
Tremaine S., 1996, AJ, 112, 105 

Henon M., 1959, Ann. d'Astrophys., 22, 126 

Henon M., 1960, Ann. d'Astrophys., 23, 474 

Henon M., 1973, A&A, 24, 229 

Hernquist L., 1990, ApJ, 356, 359 

Hernquist L., 1993, ApJ, 409, 548 

HoUey-Bockelmann K., Mihos J. C., Sigurdsson S., Hernquist 
L., 2001, ApJ, 549, 862 

Jaffc W., 1983, MNRAS, 202, 995 

Kroriawittcr A., Saglia R. P., Gerhard O., Bender R., 2000, 
A&AS, 144, 53 

Lauer T. R., Ajhar E. A., Byun Y., Dressier A., Faber S. M., 

1995, AJ, 110, 2622 
Magorrian J., Ballantyne D., 2001, MNRAS, 322, 702 
Mazure A., Capelato H. V., 2002, A&A, 383, 384 
Merritt D., 1985, AJ, 90, 1027 
Merritt D., Ferrarese L., 2001, MNRAS, 320, L30 
Meza A., Zamorano N., 1997, ApJ, 490, 136 
Osipkov L. P., 1979, Sov. Astron. Lett., 5, 42 
Palmer P. L., Papaloizou J., 1987, MNRAS, 224, 1043 
Plummer H. G., 1911, MNRAS, 71, 460 
Richstone D. O., Tremaine S., 1984, ApJ, 286, 27 
Sersic J. L., 1968, Atlas de Galaxias Australcs, Cordoba: 

Observatorio Astronomico 
Tremaine S., Richstone D. O., Byun Y.-L, Dressier A., Faber 

S. M., Grillmair G., Kormendy J., Lauer T. R., 1994, AJ, 

107, 634 

Tremblay B., Merritt D., 1995, AJ, 110, 1039 



